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We study a universal structure of two and three dimensional self-gravitating systems in the quasi¬ 
equilibrium state. It is shown numerically that the two dimensional self-gravitating system in the 
quasi-equilibrium state has the same kind of density profile as the three dimensional one. We develop 
a phenomenological model to describe this universal structure by using a special Langevin equation 
with a distinctive random noise to self-gravitating systems. We find that the density profile derived 
theoretically is consistent well with results of observations and simulations. 


INTRODUCTION 


Systems with long range forces exhibit various specific 
properties which systems with short range forces do not 
have. The prime example of them is the presence of 
another stable state differing from the thermal equilib¬ 
rium state. In this paper, we call it quasi-equilibrium 
state (QES). It has been found numerically that the dis¬ 
tribution of QES depends strongly on initial conditions 
whereas the thermal equilibrium is independent on ini¬ 
tial states. QES of Hamiltonian mean field model which 
is a toy model of systems with long range interactions 
varies from a ferromagnetic state to a paramagnetic (ho¬ 
mogeneous) state with increase of the initial total energy 
[jj. We have obtained numerically the result that the 
number density N of three dimensional self-gravitating 
system (3DSGS) in the real space can be approximated 
by the following representation when the initial virial ra¬ 
tio V(= 2K/ft) is 0 where K and H are the total kinetic 
and gravitational energy respectively Hi: 


N(r) 


(1 + r 2 /a 2 ) K 


(1) 


with n ~ 3/2, where r means the distance from the cen¬ 
ter of system. The number density of globular clusters 
consisting of about hundreds of thousands of stars which 
is the best example of 3DSGS also can be depicted by 
Eq. O with k ~ 3/2 fl (i : This density profile is uni¬ 
versal for 3DSGS. 

The recent observations have made it clear that the 
universality is not just limited to 3DSGS. As will be dis¬ 
cussed later, the cylindrical symmetric filamentary struc¬ 
ture of molecular clouds can be treated as a two dimen¬ 
sional self-gravitating system (2DSGS). The Herschel 
Space Observatory revealed the 27 filamentary structures 
in IC 5146 which is a reflection nebula in the Cygnus and 
the number density of molecular clouds have a cylindri¬ 
cal symmetry around the axis of filament Q. All of the 
number densities of molecular clouds around the axis of 
filament can be fitted by Eq. with k from 0.75 to 1.25 
jj], where r means the distance from the axis of filament. 
On the other hand, the equation m with k = 1 was uti¬ 


lized in order to describe the number density of a filament 
in the Taurus Q. So, 2DSGS has a similar universality 
depicted by Eq. m with k ~ 1. The interaction poten¬ 
tial energies of 2D and 3DSGS have a difference between 
being bounded and unbounded which will be explained 
later. However, the universality is beyond the difference. 

Obviously, the universal number density cannot be de¬ 
rived by assuming the system is in the isothermal equi¬ 
librium: The equilibrium state for 2DSGS has an exact 
solution for number density represented by Eq. m with 
k = 2 @, H). On the other hand, the number density 
of 3DSGS in the thermal equilibrium decays with r~ 2 at 
large radius r M- Therefore, a new theory explaining 
the physical mechanism behind the universality is neces¬ 
sary. 


In this paper, we shall derive the universal number 
density Eq. m from one phenomenological model par¬ 
ticular for gravity by utilizing a special Langevin equa¬ 
tion for QES and corresponding Fokker-Planck equation 
in /r dimensional space where g is 2 or 3. We shall 
treat them in the over-damped limit. Indeed, the nor¬ 
mal Langevin equation in the limit and the correspond¬ 
ing Fokker-Planck equation (i.e., Smoluchowski equa¬ 
tion) is appropriate for describing 3DSGS enclosed in 
a box near the isothermal equilibrium JT|, since the 
Maxwell-Boltzmann distribution is stable for 3DSGS in 
the thermal equilibrium state with the total mass and 
energy fixed if the radius of system is less than a criti¬ 
cal value Til, 131. However, it is well-known empirically 
that the SGSs without boundary are trapped by another 
stable state, that is QES, so that we must modify the 
Langevin and the Fokker-Planck equation in order to de¬ 
scribe the state. In particular, we shall introduce another 
noise whose foundation is so simple. 


This paper is organized as follows. Firstly, we review 
the two dimensional gravity and show that the cylindri- 
cally symmetric 3DSGS is equivalent to 2DSGS mathe¬ 
matically . Moreover, we exhibit results of N-body simu¬ 
lations of it which are the same as observations of molec¬ 
ular clouds. Next, a Fokker-Planck equation to describe 
the universal number density is derived. In doing so, we 
start a Langevin equation specialized for gravity. Lastly, 
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numerical solutions of the Fokker-Planck equation are 
shown and we reveal that the results fit well with obser¬ 
vations and numerical simulations. 


TWO DIMENSIONAL GRAVITY AND N-BODY 
SIMULATIONS OF TWO DIMENSIONAL SGS 


Two dimensional gravity 


The two dimensional gravitational potential per unit 
mass </> 2 (r) generated by mass source p 2 (r) satisfies the 
following Poisson equation in the two dimensional space. 


1 9 | d(j) 2 {r) 
r dr \ dr 


1 d 2 <f> 2 (r) 
F 2 dO 2 


47 tG' p 2 (r) 


( 2 ) 


where the subscript 2 means the two dimensional space 
and we denote the two dimensional gravitational con¬ 
stant by G' in order to distinguish it from the ordinary 
gravitational constant G. If the distribution of mass is 
circular symmetric, the above equation becomes 


1 d f dcf) 2 (r) 

-< r- 

r dr ( dr 


dirG' p 2 {r) . 


( 3 ) 


When the mass source is a mass point with a mass 
to existing in the origin, p 2 (r) = m6(0)/irr. Then, the 
Poisson equation is altered to 


A polytrope solution with a polytrope index n was 
adopted as the initial condition. The solution with 
n = oo corresponds to the thermal equilibrium state, 
so the finite n represents the deviation from equilibrium. 
Generally, the solution in the three dimensional space is 
well-known [§]. Here, we have extended it to the two 
dimensional space: The distribution function f(r,v) in 
phase space can be shown as 


f{r,v) oc 


£ n ~ x (£< 0 ) 
0 (£ > 0 ) 


( 7 ) 


where £ = (j) 2 (R) — v 2 + </> 2 (r)} and R means a radius 
of the system. Then, we run the N-body simulation by 
varying an initial virial ratio for 2DSGS V' which is dif¬ 
ferent from a virial ratio for 3DSGS (See Refs. 14] and 
0). Note that the system is not enclosed into a circle. 

Specifically when V' = 0, it is found that the number 
density in QES has a universality and the density around 
the center of the system can be fitted well by Eq. ([1]). The 
results are shown in Fig. [TJ It can be seen in Fig. QJa) 
that the number density from n = 1 to 10 has the same 
profile. In addition, figure [IJb) denotes that n which 
is the index in Eq. m ranges from 0.8 to 1.1. These 
results have a good consistency with the observations of 
molecular clouds 


FOKKER-PLANCK MODEL 


1 d f di/> 2 (r)'l 4mG'5(0) 

r dr ( dr J r 

Therefore, the potential can be solved like 

(j> 2 (r) = 2G'm\nr + const. 


( 4 ) 

( 5 ) 


Then, the interaction potential of 2DSGS is bounded 
whereas one of 3DSGS is proportional to — 1/r. 

Finally, we shall provide a brief explanation for the 
reason that cylindrically symmetric filaments of molecu¬ 
lar clouds can be regarded as 2DSGS: By using the mass 
density of molecular clouds p mc (r) and the gravitational 
potential per unit mass </> mc (r) where r means a distance 
from the axis, the Poisson equation becomes 


1 d f d</> mc (r) 

-< r - 

r dr dr 


47rGp mc (r) , 


( 6 ) 


which is mathematically equivalent to the Poisson equa¬ 
tion for 2DSGS Eq. (O. 


N-body simulations 

Because the equivalence is merely formal, we followed 
the time evolution of 10 4 -body system in the two dimen¬ 
sional space interacting by Eq. (|5J) numerically in order 
to investigate QES of 2DSGS. 


We have made sure of a fact that the number density 
of 2DSGS without boundary in QES can be represented 
well by Eq. m with k ~ 1. On the other hand, the 
number density of most globular clusters is well-known 
to be fitted by Eq. JTJ with k ~ 3/2 [l-Q . Furthermore, 
we have reported that the same density profile is obtained 
through N-body simulations of 3DSGS [i, @]. Therefore, 
we can conclude that the SGSs without boundary have 
the universal density profile depicted by Eq. ([1]) in QES. 

Here, let us derive the density profile uniformly by 
a special Fokker-Planck equation in the p dimensional 
space where p. is 2 or 3 [22j]. Before constructing the 
Fokker-Planck equation, we shall model forces influenc¬ 
ing an element of the system. That is, we shall begin by 
constructing a Langevin equation. 

We assume that the frictional force —m'yr^{t) and the 
random noise with constant intensity y/2D^(t) which 
are essential for a many-body system to reach the ther¬ 
mal equilibrium state act on the element 0 , where m is 
a mass of the element, is a position vector, and t ) 
means Gaussian-white noises |23j . The index represents 
that the vector is in the p dimensional space. Postulating 
the system to be circular (p = 2) or spherical symmetric 
(p = 3), one can see that the element is also influenced 
by a mean gravitational force —F^(r) along the radial 
direction of the system which is derived by differentiat¬ 
ing m</> M (r): — F^(r) = —md r <j^ l {r) where (f>^(r) is the 
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FIG. 1: (color online) (a) Number densities at QES for sev¬ 
eral initial polytrope indexes are plotted by open-circles. The 
curve passing the circles is a fitting function, Eq. CO- So as 
to be seen easily, each density is shifted by two digits to the 
below, (b) The optimum values of k in Eq. JlJ for fitting the 
number densities vs the initial polytope index n. 


r 




mean gravitational potential in /i dimensional space and 
r = Ir^l in which the index /i are abbreviated for sim¬ 
plify. However, this is just a mean gravity. It is natural 
to consider that the element actually is influenced by a 
fluctuating gravity around the mean value: The num¬ 
ber density producing </P(r) through the Poisson equa¬ 
tion is the mean value, and the actual distribution of 
elements must fluctuate around the value. This means 
that another noise which prevents the system from reach¬ 
ing the thermal equilibrium state is added to the normal 
Langevin equation, and so this system goes to another 
stable state, that is QES. Therefore, we can consider the 
noise distinctive to SGS. 

If assuming the intensity of the noise to be constant, 
we can obtain the following Langevin equation: 

= -F»(r){l + V2^v(t)}e? + V2D£ ,i (t) ( 8 ) 

where eis a unit vector along the radial direction in 
/r dimensional space. In the over-damped limit, equa¬ 


tion © becomes 


myr^{t) = —F^ir) (l + \/2e?y(f) j e£ 
d e 


dr 2 my 


F M (r) 2 e{f + V2D^{t) . 


(9) 


The second term on the right hand side of the above 
equation is a correction term in order to regard products 
as the Storatonovicli product 18]D 

The corresponding Fokker-Planck equation is given by 


at {my)- 


dr 2 


H — ld 
r dr 


P»(r,t) 


d—— - rP 71 {r)P^ (r,t) 
my r dr 


+ 


{my)' 2 


cf_ 

dr 2 


y-1 d , n ~ 1 d ^_ 2 


dr rf 1 1 dr 


x P M (r) 2 P M (r, t) . 


( 10 ) 


We are treating the system as a circular or a spherical 
symmetric one including N elements. So, PDF P M is a 
function of the distance from the origin r. Note that a 
relation among PDF P M , the number density AP and the 
mass density /P is as follows: P M = Af^/N = /P /{mN). 

Let us use P M (r, t) = J^{r)P^(r 7 t) instead of P 77 , 
where means Jacobian determinant, J 2 = 2nr and 
J 3 = 47T?' 2 : It can be depicting by J 77 = 2 7l_1 7rr 7 ‘ _1 . In 
doing so, we can obtain 


d D 

V^{r,t) = - - 

dt {my) z 


dr 2 


d ii- 1 


dr r 
e d 2 


V»(r,t) 


+ —^{r)V^{r,t) + 

7777 Or {m'y) z or z 


P M (r) 2 P M (r, t) . 


( 11 ) 


When the system reaches QES, d t Vq e (r) = 0. Here, 
we integrate the Fokker-Planck equation by r. Owing to 
use Pq e , the integration becomes easier: 


D 


{my) 2 
D 


+ 


eF^{r) 2 

{my) 2 

li-l 


_{my) 2 
= const. 


Ke( r ) 

2 eF^{r)F^{r) 
{my) 2 


P M (r) 


7777 


Ke(r) 

( 12 ) 


Since P qe (0) and Pq e (0) are bounded, Pq e (0) = Pq e '(0) = 
0. Therefore, the constant of the right hand side of 
Eq. (fT2l) becomes 0: 


M , f ) rF»{r) {my + 2 eF^’jr)} - (p - 1)D 

) r{D + eF»{r) 2 } qeV ’ ' 

( 13 ) 


If using the number density in QES Nq e (= NV^ e /J^), 
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we can obtain 


A^/(r) = 


rF^ir) {my + 2eF M, (?’)} + (p — 1 )eF At (r) s 
r {D + eF^(r) 2 } 


■^ e (r) 


(14) 


By substituting the gravitational potential per unit 
mass 0 A ‘(r)(= A f d rF^(r)) into the Poisson equation 
A (f>^ = = 47rt/ M 7?rA/’q e (r), an equation govern¬ 

ing F M can be obtained as follows: 

F^'(r) + ——i.F M (r) = Ai:Q^rn 2 M^ ie (r) , (15) 

where p£ e = mM^ e = mNP^ e = mNV^ e / J M and 
relates quantities appeared previously as Q 2 = G' and 
G 3 = G. 

Here, we nondimensionalize these equations by using 
the following units of length and force: 


[length] = J 


_ m 2 (m+2 )d _ 

27rQF-m 2 J\fq e (0){i-L 2 m / ‘y-\-8n(fj, 2 -\-4fi-\-2)eQF-m 2 Afq e (0)} 

(16) 


[force] = 


871-A(^+2)-Pg M ^ 2 ACe(0) 

// 2 //io +S7t(// 2 -(-4//+2yC?''m 2 A ; 7A f :) 


Then, equations (fill) and (fl5l) are altered to 


(17) 


Ke(r) = 

fF»(f){l + 2pqF^(f)}+p(p-l)qF»(f) 2 
^ > f{p + 2(p 2 +4p + 2)q + 2p 2 (p + 2)qF»(f) 2 } 

( lg ) 


and 

^'(f) + ^^(f)=^ e (f), (19) 

where q = 47reG /J 'm 2 Afq e (0)/(pmj ) and variables with 
over bars denote dimensionless. We should solve these 
equations with boundary condition Afq e (0) = 1. 


RESULTS 


although the deviation between two curves increases as 
r gets larger. Therefore, we can understand that the 
best-fit curves are derived from our model by varying q 
appropriately. 




FIG. 2: (color online) Numerical solutions of Eqs. (1181) and 
Gi) for p = 2 (a) and p = 3 (b). As the curve changes from 
the left to the right in (a), q gets larger from 0 to 0.28 in steps 
of 0.07. On the other hand, q gets larger from 0 to 0.05 in 
steps of 0.01 in (b). The (red) dashed curve in (a) and (b) 
means 1/(1 + F 2 ) and 1/(1 + r 2 ) 3 / 2 , respectively. 


Finally, we examine the range of the index n. From 
Eqs. (fl8l) and (fliil) . the number density around the center 
of the system in QES can be described by 


ACM 


1 

(1 q. f2)«( 9 ) 


( 20 ) 


The numerical solutions for p = 2 and 3 are shown in 
Fig. [2] by changing q. The curves with q = 0 on both 
figures correspond to the thermal equilibrium state. For 
comparison with observations, l/(l+r 2 ) and l/(l+r 2 ) 3 / 2 
are also plotted by a dashed curve in Fig. [2] (a) and 
(b), respectively. The dashed curves are typical best-fit 
ones for densities of molecular clouds or globular clus¬ 
ters. From Fig. [2] (b), one can see that the numerical re¬ 
sult of our model with q = 0.01 completely coincide with 
the typical number density. Fig. [2] (a) also shows the 
good agreement of our model with observations and nu¬ 
merical simulations for small radius by setting q = 0.56, 


where n is a function of q: 

/ \ _ + 2 ) {1 + + l)g} 

p + 2{p 2 + 4 p + 2 )q 


( 21 ) 


From this equation, because of q > 0, one can see the 
range of k as follows: 


(p + 1 ){p + 2) < K< V + 2. 

2 (p 2 + 4 p + 2) ~ p 


( 22 ) 


Therefore, with regard to p = 2 and p = 3, | < k < 2 
and || < k < §, respectively. Both ranges include the 
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observed indexes. However, the observed ranges of the 
index are so narrower, which means that the value of q 
is limited. The limited q can be regarded as a new kind 
of fluctuation-dissipation relation [19j which is particular 
for SGS because q includes a ratio of the intensity of 
mean gravity fluctuation e to the friction coefficient 7 . 
Indeed, this phenomenological approach cannot provide 
an answer to a question why the specific value of q is 
selected. 


CONCLUDING REMARKS 


In this paper, we made it clear numerically that 2DSGS 
without boundary goes to QES in which the density pro¬ 
file is depicted by Eq. © with k ~ 1. It is well-known 
that QES of 3DSGS without boundary can be described 
by the same one through the observations of globular 
clusters and N-body simulations. That is, it was shown 
that there is the universal QES for the SGSs: As men¬ 
tioned before, the interaction potential at 2DSGS are 
qualitatively different from one at 3DSGS. Moreover, the 
isothermal equilibrium of SGS with boundary differs be¬ 
tween 2D and 3D loj, 11, 2fl 211. On the contrary, the 
paper unveiled that 2D and 3DSGS without boundary 
have the same universality in QES. Furthermore, we de¬ 
veloped the model to derive the universal density pro¬ 
file from the special Langevin equation including the dis¬ 
tinctive noise of SGS. Indeed, the solution of the corre¬ 
sponding Fokker-Planck equation in QES was depicted 
by Eq. ©• In addition, it was found that the index k in 
Eq. ([1]) can be represented as a function of the intensity 
of the particular noise for SGS, the friction coefficient, 
and others. So, we showed the range of n for each di¬ 
mension analytically which has a good consistency with 
observations and simulations. 

The theory justifying the frictional force and the ran¬ 
dom noise is indispensable for constructing the Langevin 
equation. As mentioned before, there are theories for 
3DSGS [H, 17|, but not for 2DSGS. Therefore, such a 
theory must be developed. 
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